A Quantum Top Inside a Bose Josephson Junction 
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We consider an atomic quantum dot confined between two weakly-coupled Bose-Einstein conden- 
sates, where the dot serves as an additional tunneling channel. It is shown that the thus-embedded 
atomic quantum dot is a pseudospin subject to an external torque, and therefore equivalent to a 
quantum top. We demonstrate by numerical analysis of the time-dependent coupled evolution equa- 
tions that this microscopic quantum top is very sensitive to any deviation from linear oscillatory 
behavior of the condensates. For sufficiently strong dot-condensate coupling, the atomic quantum 
dot can induce or modify the tunneling between the macroscopic condensates in the two wells. 
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In the field of ultracold atoms, whose most spectacu- 
lar achievement on relatively large scales is Bose-Einstein 
condensation (BEC), not only macroscopic systems are 
of interest, but also to confine several or even single 
atoms into optically created microtraps is becoming a po- 
tentially important experimental tool of what might be 
coined "nanobosonics." In nanoelectronics the control of 
electronic quantum dots is performed by biased conduct- 
ing leads, attached to it. In nanobosonics the role of the 
"leads" is played by finite superfluid reservoirs of given 
particle number, which can be coupled to a particular 
atom by optical transitions. Trapping and manipulating 
single atoms [J, Q] opens up new perspectives in the co- 
herent control of quantum states, and is of relevance for 
quantum computational tasks 

It has recently been demonstrated by Recati et al. [4( 
that an atomic quantum dot (AQD) (a single atomic 
two-level system), optically coupled to a superfluid BEC 
bath, can be mapped onto the spin-boson model. This 
system then exhibits a dissipative quantum phase transi- 
tion, characteristic of this model 0|. Here, we study such 
a spin-boson model, but with a time- dependent bath: An 
AQD located inside a Bose Josephson junction (BJJ), 
i.e., a single bosonic atom coupled to two superfluid reser- 
voirs. The setup under consideration is schematically de- 
picted in Fig.[TJ The Bose-Einstein condensate is trapped 
by the double- well potential Vbec(t)- The atom of the 
dot, which is in a hyperfine state different from that of the 
condensate, is confined by a very tight potential Vaqd, to 
which condensate atoms are insensitive, and which causes 
a large gap for double occupation of the dot. The cou- 
pling of the dot to the condensates in the wells is per- 
formed in a tunable way via a Raman transition Q. Due 
to their coherent nature, the weakly-coupled condensates 
exhibit quantum tunneling [fj]. In the present paper, we 
investigate the mutual influence of the induced conven- 
tional Josephson oscillations between the two wells and 
the AQD, which provides an additional tunneling chan- 
nel. We demonstrate, by numerically solving the time- 
dependent coupled evolution equations of AQD and con- 
densates, that this additional channel can in certain cases 



directly affect the macroscopic Josephson tunneling. 
The Hamiltonian of our system consists of three parts 



H — -ffcond + -ffdot + -ffcoupl- 



(1) 



We will first describe these three parts in turn. The part 
-Hcond characterizes the double- well trapped BEC: 



J dr |V(r,i) 



2m 



v 2 + y BEC (r) 
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&(r,t) 
(2) 



where &(r,t) is the condensate wavefunction and m the 
atomic mass. We assume that at low energies the inter- 
particle interaction is given by the usual pseudopotential 
V(r — r') = g5(r — r'), where g — Ai:h 2 a s /m, and a s 
is the s-wave scattering length. The condensate is de- 
scribed within Gross-Pitaevski! theory, sufficiently accu- 
rate at the very low temperatures we are considering 0]. 

For the present dilute bosonic gas of finite extent, the 
quantum tunneling between the two wells is adequately 
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FIG. 1: [Color online] An atomic quantum dot between two 
weakly-coupled condensates, trapped in a double-well poten- 
tial Vbec- The dot is a simple two-level system = single 
atom present/not present and is created by the tight potential 
Vaqd, located at the position of the top of the barrier. Atoms 
can be exchanged between wells either by direct tunneling 
(dashed arrows) or via the dot, coupled to the condensates by 
a transfer matrix T. 
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described within a two- mode approximation 0, 0|; one 
expands V BEC (r) around each minimum, and introduces 
the local mode solution, (j>x i2 (t), for each well sepa- 
rately. In first approximation the two modes can be 
considered to be orthogonal, J drcj)i(r)(f) 2 (r) — 0. The 
two-mode approximation then results in the following 
ansatz for the total condensate wavefunction (3^,2 (t) = 



(3) 



Since we are interested in tunneling events, i.e., in the 
time dependence of the wavefunctions \P\ j2 (t), it is con- 
venient to write -ff C ond in the form 



ffcond - J2 E°\m)\ 2 + u t \m)\ 



i=l,2 



-K^(t)^(t)+^(t)^(jt)), (4) 



dr 



where E? = f [-£ | Vfc(r) | 2 +\Ur)\ 2 V BEC (r) 
are the zero-point energies in the wells 1 and 2, re- 
spectively, the effective "on-site" interaction between the 
particles is given by Ui = g J \<fii(r)\ 4 dr, and finally 

K = -Sl£&Mr)V<h(r))+Mr)V BEC (r)Mr)]dr de- 



(5) 



notes the coupling matrix element 

The Hamiltonian of the dot itself is given by 

ffdot - Jdr[-h5$(r,t)d{r,t) 



U dd 



S{r,t)d\r,t)d(r,t)d(r,t)}. 



We assume that the dot operator factorizes according 
to d(r,t) — d{t)(f>d{r), where 4> d {r) is the spatial wave 
function of the atom on the dot normalized to unity, 
J dr\(f) d (r)\ 2 = 1. The repulsive interaction between the 
dot atoms we consider to be much larger than any other 
energy scales in the system, Udd — * oo. The dot can then 
be described as a two-state system, the two states being 
that an atom is or is not trapped inside the dot. Finally, 
the dot interacts with the condensate as follows 



H coup i=g d c / dr\<P(r,t)\ 2 S(r,t)d(r,t) 



+ hfl dr(&*(r, t)d(r, t) + h.c). 



(6) 



Here, g dc is the dot-condensate interaction constant, and 
the second term describes the coupling of the condensate 
atoms to the lowest vibrational state in the AQD via a 
Raman transition with characteristic Rabi frequency fl. 
Spontaneous emission is suppressed by a large detuning 
from the excited electronic states, which is absorbed into 
the effective dot energy K5 

To represent the evolution equations following 
from the Hamiltonian P in a physically trans- 
parent form, we introduce a new set of pa- 
rameters U id = gdc J dr\4>i(r)\ 2 \4> d {r)\ 2 , U 12d = 



gdc J dr0i(r)0 2 (r)|(/)rf(r)| 2 , and T, = Ml f dr(f>i{r)(j> d (r), 
with Ui2d — U* 2d and Tj = T* . In the single-occupation 
limit, the temporal wavefunction of the dot is just a su- 
perposition of singly and non-occupied states, \&d(t)) = 



ao(t)\0) where |ao| 



|ai| 2 = 1. The dot op- 



erators then correspond to Pauli matrices: d(t) — > <r_ (t) 
and <v(t) — > a+(t), introducing the pseudo-spin ladder 
operators a± = \{<J X ± i& v ) in terms of the Pauli matri- 



ces cr. 



We can thus write for the coupling term 



Hcoupl = J2 [ U id\*i\ 2 + (Ul2d#!*2 + h, 
i=l,2 

+ {r i *jo-+(t) + h.c} 



(7) 



One can now derive the coupled equations of motion for 
the condensate © and the spin s(t) = (W d (t)\&\W d (t)) = 
(W d \a{t)\W d ), from the total Hamiltonian Eq. |(T]). The 
equations for the condensate are 



ihdt*i = + U x Ni(t) + U ld n d (t)] <fi 
+(Ui 2d n d (t) - k)&2+T!S-, 

ihd t V 2 = [El + U 2 N 2 (t) + U 2d n d (t)] $ 2 
+ (U 12d n d {t) - k)&i+T 2 S-, 



(«) 



while the dot equations are 



ifidts- = {-h5 + U u Ni(t) + U 2d N 2 {t) + U 12d &*& 2 

+U 12d n^i\ s- - {Ti&x + T 2 V 2 )s z , 
ihd t s z = 2(Ti«f + T 2 <? 2 *)s_ - 2(T 1 !Pi + T 2 W 2 )s+. (9) 



It is easily verified that the Eqs. ( 
vector form of a Bloch equation 



can be written in the 



hdts = u>(t) x s, 
where the time-dependent frequency vector reads 



(10) 



2Ti y/ mjt) c os Bi (t) + 2T 2 y/ N 2 (t) c os 6 2 (t) 
ujH)= j -2TiV^iWsin0i(*) -2T2V / ^2Wsin0 2 (i) 
-HS + UuN^t) + U 2d N 2 (t) + u 12 (t) cos<p(t) 

, (11) 

where to 12 {t) = 2U 12d y/N 1 (t)N 2 (t) and 4>{t) = 6 2 {t) - 
6i(t). It follows that the AQD inside the BJJ is equiva- 
lent to a quantum top. In the case of time-independent 
uj, Eq. iflOl) can be solved analytically. The presence of 
Josephson tunneling between the condensates however 
generally results in a time-dependent <*> = u>(t), and the 
equations need to be solved numerically [id ]. 

For simplicity, in what follows we consider the case of 
a fully symmetric system: E± — E 2 = 0, U\ = U 2 = U , 
U\d = U 2d , Ti = T 2 = T, Ui 2d = U 2 i d . In order to com- 
pare our results with previous work on BJJ |9], we intro- 
duce dimensionless parameters: t — » 2nt, A = UNq/k, 
where iVo = A^i(0) + ^(0) is the initial total number 
of particles in the condensates; note that the conserved 
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population imbalance 

»(*) = 



N 1 (t)-N 2 (t) 

Nn 



(12) 



FIG. 2: Results for weak coupling T rc i = 0.01 (N = 1500 
throughout Figs. 2-4, as used in experiment [l3y) and A = 10. 
n(0) = (Fig.2a - solid line), n(0) = 0.5 (Fig. 2a - dashed 
line). Fig. 2b displays the dot occupation for n(0) = 0; Fig. 2c 
shows the precessional behaviour of u> (in units of 2k), and 
Fig. 2d the corresponding spin nutation for n(0) = 0. 




FIG. 3: Condensate in the self-trapped MST state, n(0) = 
0.8, A = 10, T ro i = 0.01. Relative population oscillations are 
shown in Fig. 3a, and the dot occupancy in Fig. 3b. In Fig. 3c, 
we display the projection of the frequency u> lfT0|) on the x — y 
plane in units of 2k, and in Fig. 3d the pseudospin. 



quantity is N to t — Ni(t) + N 2 {t) + na{t). For numeri- 
cal convenience, we fix the dot energy at M = 2k. In 
addition, the interactions between AQD and condensate 
Uid/i^n) and ?7i2d/(2/c) are assumed to be vanishingly 
small (Uid <C U, Uud <C U), and lu z ~ — M = const. Our 
main parameters are then the dimensionless strength of 
coupling of dot to condensate T ro i — T/k, quantifying the 
relative importance of tunneling channels via dot and di- 
rectly by conventional Josephson tunneling, respectively; 
and A, measuring the relative importance of mean-field 
interaction in and tunneling between the wells. 

In the following, we present results for the fractional 



the occupation of the dot n<i{t) and the trajectories of 
the pseudospin s on the Bloch sphere, as well as the 
projection of the frequency-vector uj on the x — y plane. 

We first consider the situation when the dot does not 
have a notable effect on the tunneling between the wells 
(Figs. 2 and 3); we fix A = 10, and only change the initial 
condition for the particle imbalance, n(0). The most sim- 
ple situation is the stationary one of an initial population 
imbalance n(Q) = and initial phase difference (f)(0) = 
(for definiteness in all figures n<j(0) = 1, i.e., there is 
initially exactly one atom in the dot). These conditions 
result in an AQD coupled to a time-independent BEC 
i.e., to the problem of a spin in a constant magnetic field, 
however without dissipation. The pseudospin generally 
undergoes nutation (also if we put n(0) = 0.5 - Fig. 2a 
- dashed line), as shown in Fig. 2d, while the vector u> 
precesses, Fig. 2c. However, there is an exception to this 
general behavior: For A = 1, n(0) = there occurs a sim- 
ple precession of the pseudospin (not shown), while 
the occupation of the dot exhibits linear oscillations. For 
A / 1 the precession is lost, an effect due to the finite 
number of particles in the system. 

The fact that Ntot is a finite quantity constitutes one 
major difference to the system of a singl e sp in coupled 



to superconducting leads considered in [11]. Further- 
more, while deviation from simple precessional behavior 
also occurs in that system, the effective fractional pop- 
ulation imbalance n(t) is essentially zero. Regimes re- 
lated to large n(t) of order unity, to be discussed below, 
are thus not accessible for the superconducting Joseph- 
son junction - single spin system. In addition, the tun- 
neling (quasi-)particles are treated as noninteracting in 
the latter case. Here, by contrast, including interactions 
between the fundamental bosons is crucial. In particu- 
lar, as a consequence of interactions, and as discussed 
in detail in @, E|, depending on A and the initial con- 
ditions, a condensate in a double-well potential can ex- 
hibit a novel quantum state - macroscopic self-trapping 
(MST), successfully observed experimentally fUj. MST 
is only present for the self-interacting matter waves, and 
is characterized by a nonzero time average of the popu- 
lation imbalance n(t). The transition to the MST state 
is a gradual crossover, and we observe that our quantum 
top is very sensitive to this crossover. In the plots of 
Fig. 2, far away from the self-trapped state, the coupling 
strength T re \ does not influence in a qualitative way the 
behavior of the quantum top. The pseudospin behavior 
however drastically changes as we approach MST It ap- 
pears that the AQD is sensitive to the deviation from 
linear oscillatory behavior of the condensates occurring 
in this regime. The linear oscillation of the dot occu- 
pation is then destroyed (Fig. 3b), and the pseudospin 
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FIG. 4: [Color online] Results for a 7r-junction. In Fig. 4a, 
self-trapped MST state with ra(0) = 0.6, A = 1.25; decoupled 
dot, T IB \ = (black dashed line); T rc i = 0.1 (thick red line), 
Trci = 1 (thin wavy blue line); the pseudospin then nutates, 
Fig. 4c (T r ei = 0.1). The AQD induces strong modifications 
both for small oscillations between wells and in the MST state, 
Figs. 4b and 4d; n(0) = 0.01, and T Icl = (dotted curve), 
Trei = 1 (solid curve). Weak interaction coupling, A = 0.1 in 
Fig. 4b, strongly coupled MST state, A = 1.1, in Fig. 4d. 



undergoes multiple-frequency rotations (Figs. 3c and 3d). 
In the MST state, the pseudospin can thus behave in a 
rather irregular manner already for small values of the 
relative coupling T rc \. 

There is another potentially interesting regime, which 
occurs when the effect of the dot on the tunneling be- 
tween the wells becomes significant (Fig. 4), i.e., with 
increasing value of the couplingto the wells T re i. Con- 
sider, for instance, a 7r-junction [9J, <fi(0) = it, n(0) = 0.6, 
A = 1.25 (Fig. 4a, dashed line). The coupling to the dot 
leads to small oscillations between the wells (results for 
different T ro i are shown in Fig. 4a), and the pseudospin 
undergoes nutation, as apparent from Fig.4c. For weak 
interactions (in the so-called Rabi regime and very 

small particle imbalance, the effect of the dot becomes 
more pronounced (Fig. 4b) . When the coupling to the dot 
is very weak, we observe the nutation of s and precession 
of UJ. Increasing T re i leads to significant modifications of 
the tunneling picture (Fig. 4b, solid line), with strongly 
non-sinusoidal oscillations of the population imbalance. 
Finally, we observe that, changing T rc \ from small to large 
values, the dot can switch the BJJ from the MST state 
to a small population imbalance state (Fig. 4d). 

In conclusion, we have shown that two weakly-coupled 
condensates, with an AQD situated at the location of 
the top of the barrier between them, can exhibit several 
regimes of oscillatory behavior. The AQD behaves as a 
quantum top whose behavior is very sensitive to the tun- 
neling mode between the condensates. Even for small 
couplings and stationary condensates the "spin" of the 
dot nutates, an effect due to the finite number of parti- 



cles in the system, which vanishes for an infinite system. 
Nutation is a characteristic feature of the quantum top in 
regimes far away from the MST state. Conversely, mov- 
ing towards the self-trapped regime, we obtain strong de- 
viations from nutational behavior, and the quantum top 
motion becomes strongly irregular. However, when the 
AQD itself modifies in a significant way the oscillations 
between the wells, nutation can emerge also in a MST 
state. Finally, the dot can act as a switch for the BJJ 
from MST to small population imbalance oscillations. 

We treated the condensate on a mean-field level. In 
future studies, it would be of interest to study the in- 
fluence of condensate quantum fluctuations on the AQD 
in the limiting case that the dot provides the domi- 
nant tunneling channel between the condensates. 
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